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Abstract: We perform a lattice analysis of the Faddeev-Niemi effective action conjectured 
to describe the low-energy sector of SU{2) Yang-Mills theory. To this end we generate 
an ensemble of unit vector fields ('color spins') n from the Wilson action. The ensemble 



does not show long-range order but exhibits a mass gap of the order of 1 GeV. From 
the distribution of color spins we reconstruct approximate effective actions by means of 
Qh! exact lattice Schwinger-Dyson and Ward identities ('inverse Monte Carlo'). We show that 

the generated ensemble cannot be recovered from a Faddeev-Niemi action, modified in a 
minimal way by adding an explicit symmetry-breaking term to avoid the appearance of 
^ " Goldslono modes. 
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1. Introduction 

Confinement in pure Yang-Mills theory is still a theoretical challenge. The problem ac- 
tually has two faces. On the one hand, there is confinement of static external sources in 
the fundamental representation which manifests itself through the appearance of a linear 
potential (nonzero string tension). On the other hand, there should be gluon confinement 
implying a finite range of the gluonic interactions, i.e. a mass gap. How the two faces are 
related is largely unclear at the moment. 

Recently, Faddeev and Niemi (FN) have suggested that the infrared dynamics of glue 
might be described by the following low-energy effective action , 

m 2 (^n) 2 + ^H^H^ . (1.1) 

Here, n is a unit vector field with values on S 2 , n 2 = n a n a = 1, o = 1,2,3; m is a 
dimensionful and e a dimensionless coupling constant. The FN 'field strength' is defined 
as 

Hfj, v = n ■ d^n x d u n . (1.2) 



Sfn 



ctx 
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Faddeev and Niemi argued that (1.1) "is the unique local and Lorentz-invariant action 
for the unit vector n which is at most quadratic in time derivatives so that it admits a 
Hamiltonian interpretation and involves all such terms that are either relevant or marginal 
in the infrared limit" . 

It has been shown that SpN supports string-like knot solitons [||, ||, |j, characterized 
by a topological charge which equals the Hopf index of the map n : S 3 — > S 2 . Here, 
n is supposed to be static and approaches a uniform limit at spatial infinity, = e z . 
In analogy with the Skyrme model, the H 2 term is needed for stabilization. The knot 
solitons can possibly be identified with closed gluonic flux tubes and are thus conjectured 
to correspond to glueballs. For a rewriting in terms of curvature-free SU(2) gauge fields 
and the corresponding reinterpretation of 5fn we refer to ||. 

In order for the model to really make sense, however, the following problems have to 
be solved. First of all, neither the interpretation of n nor its relation to Yang-Mills theory 
have been fully clarified. An analytic derivation of the FN action requires 

• an appropriate change of variables, A — > (n,X), relating the Yang-Mills potential A 
to n and some remainder X 

• the functional integration over X to arrive at an effective action S e g for the n-field. 
Some progress in this direction has been made || |7], ^, |9|, [l^, [ll| on the basis of the 



Manton-Cho decomposition |L3, 13], 



= C^n - n x d^n + W ^ , (1.3) 

where C is an Abelian connection and n ■ = 0. Nevertheless, it is fair to say that there 
are no conclusive results up to now. 

Second, there is no reason why in a low-energy effective action for the n-fields both 
operators in the FN 'Skyrme term', which can be rewritten as 

H 2 = (d^n ■ d^nf - (d^n ■ d u n) 2 , (1.4) 

should have the same coupling. Third, and conceptually most important, Sfn has the 
same spontaneous symmetry breaking pattern as the nonlinear cr-model, SU(2) — > U(l). 
Hence, it should admit two Goldstone bosons and one expects to find no mass gap. In 
order to exclude these unwanted massless modes we have recently suggested to break the 



global SU(2) explicitly [14|, an idea that has subsequently also been adopted by Faddeev 
and Niemi [jl5|] . 

In what follows the FN hypothesis will be tested on the lattice. To avoid the appearance 
of Goldstone bosons we allow for explicit symmetry-breaking terms. 

2. Generating an SU(2) lattice ensemble of n— fields 

The conceptual problem to be solved in the first place is to obtain a reasonable ensemble 
of n-fields. The (lattice version of the) decomposition ( |1.3| ) is of no help: it assumes some 
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particular choice of n on which the decomposition is then based. One way of defining an 
ra-field is via Abelian gauge fixing, originally introduced by 't Hooft pi. A prominent 
example in this class of gauges is the maximally Abelian gauge (MAG) which is obtained 
via maximizing the functional [FT?]] 



*mag[^; 9] = tr ( T 3 9U ^ T 3 9U 1^ =Y, tr {^U^n^Ul^ = F MAG [U; n] , (2.1 



X,fl x,n 

with respect to the gauge transformation g. The maximizing g then defines the n-field 
according to 

n x = gt T 3 g x = n x -r. (2.2) 
Instead of maximizing -Fmag with respect to g one can equivalently maximize .Fmag with 



respect to n [18] which results in the condition 

A[U]n x = X x n x . (2.3) 

Here, A[C7] denotes the covariant Laplacian in the adjoint representation (see App. |A[ ), 
while X x is a Lagrange multiplier imposing that n x is normalized to unity, locally at 
each lattice site x. In principle, ( [2.3| ) can be solved for the field n associated with the 
background U. However, as this background is distributed randomly along its orbit we will 
in turn obtain a random ensemble of n-fields characterized by the two-point function 

G% = «<) = U ab 5 xy . (2.4) 

Thus, nontrivial correlations are absent. Fig. |l| shows that this is indeed what one gets in a 
typical Monte Carlo run. One way out of this problem is to follow the continuum approach 
of H which starts out with a covariant gauge fixing. After having generated SU (2) lattice 
configurations using the standard Wilson action we therefore fix to lattice Landau gauge 
(LLG). The latter is defined by maximizing the functional 

i^LLG = ^2tr n U XjlM (2.5) 
with respect to the gauge transformation S7. In this way we impose some 'preconditioning' 



| IS] which (i) eliminates the randomness in our Yang-Mills ensemble and (ii) leaves a 
residual global SU (2)-symmetry. The Landau gauge configurations are then plugged into 
the MAG functional (|2.1[ ) which subsequently is maximized with respect to g. The gauge 
transformation g obtained this way determines n according to ([2.2;). One may say that 
g (and hence n) measure the gauge-invariant (!) distance between the LLG and MAG 
gauge slices (see Fig. ||). In App. [B] we show that LLG and MAG are 'close' to each other. 
Therefore, the maximizing g is on average close to unity, hence, on average, n will be 
aligned in the positive 3-direction. In this way we have explicitly broken the global SU (2) 
down to a global U(l). 

All computations have been done on a 16 4 -lattice with Wilson coupling (3 = 2.35, 
lattice spacing 0.13 fm and periodic boundary conditions. For the LLG we used Fourier 
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Figure 1: Behavior of the two-point function G^L along a lattice axis for a random ensemble of 



n-fields, obtained via MAG and (2.2). Note that also the value 1/3 for a = b, x = and x = L is 
correctly reproduced. 
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Figure 2: Gauge invariant definition of n = g^r^g. The gauge equivalent configurations A\ and 
A2 are both mapped onto the same 'representatives' on the LLG or MAG slices (ignoring Gribov 
copies). Thus, they are both associated with the same gauge transformation g defining n. 



accelerated steepest descent |2(| (see Fig. ||). The MAG was achieved using two indepen- 
dent algorithms, one (AI) using iterations based on elementary geometric manipulations 
(including overrelaxation steps), the other (All) being analogous to LLG fixing (see Fig. |I]). 
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Figure 3: Behavior of the LLG-functional using different algorithms. The parameter measures 
the 'distance' from the LLG, i.e. for = the LLG is achieved. 
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Figure 4: Behavior of the MAG-functional using different algorithms. 
3. Numerical results 

As expected, we observe a non-vanishing expectation value of the field in the 3-direction, 
a 'magnetization' 50? defined through (n a ) = DJl5 a3 . Thus, the global symmetry is in- 
deed broken explicitely according to the pattern SU(2) — > U(l). We demonstrate this by 
exhibiting the angular distribution of the n-field on its target space S 2 in Fig. |. The 
azimuthal angle <p is equally distributed, while the distribution of the polar angle 6 has a 
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maximum near ir/2 corresponding to the north pole, n = (0,0, 1). 
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Figure 5: Distribution of polar and azimuthal angles (9 and <j>) associated with the unit vector n 
on S 2 . The uniform distribution for <fi and the maximum for 9 = ir/2 shows that n is located near 
the north pole, n = (0, 0, 1). 

Explicit symmetry breaking also shows up in the behavior of the two-point functions 
(Fig. |6|). The longitudinal correlator, g\ = (n^n 3 ,) ~ (re 3 ) (re 3 ) = 9J1 2 , exhibits clustering 
for large distances, the plateau being given by the magnetization (squared). The transverse 
correlation function (of the would-be Goldstone bosons) 

1 2 

G^ = G^ = -J2(n>i), (3.1) 

i=i 

decays exponentially as shown in Fig. ^. This means that there is a nonvanishing mass gap 
M whose value can be obtained by a fit to a cosh-function (see Fig. [?]). The numerical 
values of the observables, Wl, M and the transverse susceptibility, 

X 

which all can be derived from the two-point functions, are summarized in Table 1 for 
both algorithms. The disagreement between AI and All is statistically significant. We 
attribute it to the ubiquitous Gribov problem (for Abelian gauges, see 
the lattice, this is the statement that maximizing gauge fixing functions like -Fmag or Fllg 
is equivalent to a spin-glass problem with an enormous number of degenerate extrema. 
This implies that the algorithms AI and All will almost certainly end up in different local 
maxima, which explains the difference between rows one and two in Table 1. 

As shown in the last column of Table 1, the numerical results for the mass gap M lead 
to a value of about 1 GeV in physical units. 
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Figure 6: Behavior of the two-point correlators of the n-field along a lattice axis (labelled by 
coordinate x). Note the difference between algorithms AI and All (FFT). Error bars exhibit the 
statistical error of the Monte Carlo simulation. 
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Figure 7: The transverse correlation function along an arbitrary lattice axis, fitted as G ± (x) = 
ai cosh(M(a; - L/2)) + a 2 with on = 0.0048, a 2 = -0.0053, M = 0.6084. Data points are obtained 
with algorithm AI. 



n 3 is a local functional of the n\ n 3 = (1 — n l n t ) 1 ^ 2 . Thus, one expects the same expo- 
nential decay for the longitudinal correlator G'L This can be confirmed with a numerical 
value for the mass gap of M = 0.66 a -1 . 
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algorithm 


£01 


-4 _L 


aM 


M [GeV] 


AI 


0.438 


92.57 


0.61 


0.95 


All 


0.366 


79.66 


0.67 


1.03 



Table 1: Numerical results for some observables as obtained from the longitudinal and transverse 
two-point functions, G" and G -1 , respectively. 

x 10- 3 1 1 1 1 1 1 1 1 1 
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Figure 8: The time-slice correlator fitted to a (properly symmetrized) Bessel function, C if) = 
c x [t 2 K 2 {Mt) + (t- L) 2 K 2 {M{L - t))] + c 2 , where a = 0.0006, c 2 = 0.0001 and M = 0.6423. Data 
points are obtained with algorithm AI. 

To improve statistics, we have calculated the time-slice correlator, 

C\t) = L~*Y,Git ■ (3.3) 

X 

In the continuum, for purely exponential decay of G -1 -, this would become proportional to 
a modified Bessel function K.%. An associated fit works very well as is shown in Fig. |8|. 
Fitting the time-slice correlator according to Fig. ||, we obtain for the mass gap 

aM = 0.642 i.e. M = 0.97 GeV . (3.4) 

This is the value with the smallest statistical errors. 

The mass gap obtained differs significantly from the SU{2) mass gap, M S jj^ — 1.5 
GeV, obtained directly from a Wilson ensemble with (3 = 2.4 J24[ |. We believe that the 
difference is due to the highly nonlocal relation between the original Yang-Mills degrees of 
freedom (the link variables) and the color spin n. After all, we have implicitly solved the 
partial differential equation ( p.3[) with link variables U in Landau gauge entering the adjoint 
Laplacian. The solution n will clearly be a nonlocal functional of these f/'s. Consequently, 



-8- 



we cannot expect that the exponential decay of G^fn] will be governed by the lowest 
excitation of the [/-ensemble. 

4. Effective action and Schwinger— Dyson equations 

At this point it is natural to ask whether there is an effective action Sefffn] that reproduces 
the distribution of n-fields leading to the results of the previous section. 

At low energies, it should make sense to employ an ansatz in terms of a derivative 
expansion, 

S cS = £ *jSj[n] + E X 3 S >' h \ ' ( 41 ) 

3 3 

with 0(3) invariant operators Sj and noninvariant operators Sp which are ordered by 
increasing mass dimension. Up to dimension four, one has the symmetric terms, 

S l = (n,An), S 2 = (n,A 2 n), 

S3 = (n ■ An, n • An) , S4 = (n • djtd^n, n • c^c^n) , 



and the symmetry-breaking terms including a unit vector 'source field' h |l4j (which can 
be thought of as the direction of an external magnetic field), 

S[ = (n,h), S' 2 = {n-h,n-h) , S' 3 = (n ■ An, n ■ h) . (4.3) 

In the above, we have introduced the scalar products 

(f,9) = J2f*9*> u-v = u a v\ (4.4) 

X 

and the usual lattice Laplacian A (see App. [A]). 

Note that the n-field configurations are classified by the Hopf invariant irrespective 
of the particular form of the (effective) action. This, together with the usual scaling argu- 



ments, shows that the action (4.1) with the operators (4.2) and (O) should still support 



classical knot soliton solutions. Our ansatz thus does not exclude this important feature. 



The couplings in (4.1) can be determined by inverse Monte Carlo techniques. The no- 
tion is suggestive: instead of creating an ensemble from a given action, one wants to com- 
pute a (truncated) action which gives rise to the given ensemble of n-fields. A particular 
approach uses the Schwinger-Dyson equations [25, |2q |. These represent an overdetermined 
linear system which can be used to solve for the couplings in terms of correlation functions. 
The latter are nothing but the coefficients of the linear system. 

For an unconstrained scalar field <fi, the Schwinger-Dyson equations follow from trans- 
lational invariance of the functional measure, implying 

= j V4>P x {F[4>]exp(-S[4>])} , (4.5) 

where P x = —iS/5(j) x is the (functional) momentum operator, and F an arbitrary func- 
tional 1 of the field (\>. For a constrained field like n with a curved target space things are 

Usually one chooses F[cj>] = 4>(xi) . . . <j)(xk). 
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slightly more subtle |25f| . There is, however, a rather elegant way to derive the Schwinger- 
Dyson equations if one exploits the isometries of the target space S 2 f27|. The target space 
measure, 

Vn = Yl dn x 5(n 2 x - l) , (4.6) 

X 

is obviously rotationally invariant, i.e. under n — > Rn, R £ 0(3). This implies the modified 
Schwinger-Dyson identity 

J VnL x {F[n]exp(-S^[n,h])} = 0, (4.7) 
where L x denotes the angular momentum operator (at lattice site x), 

iL x = n x x— or iVt = e abc n b x — . (4.8) 
dn x x x bn% v ; 

In shorthand-notation, ( |4.7| ) can be rewritten as 

(L x F[n) - F[n)L x S eS [n]) = . (4.9) 



These exact identities can be used to determine the unknown couplings A-,-. To this end 
one chooses a set of field monomials Fi [n] and plugs them into (^1]) together with the form 
Q4.1| ) of the action. This yields the local linear system 

J2(Fi L x Sj)\j + L x S' k )X' k = (L x Fi) , (4.10) 

j k 

which, in principle, can be solved numerically, for instance by least-square methods. The 
identities obtained so far hold for arbitrary actions S^n]. In particular, we have not made 
use of any symmetries. Taking the latter into account will lead to Ward identities. 

Let us specialize to our lattice effective action (4.1). It is a sum of a symmetric part 



S containing the terms ( |4.2[ ) and an asymmetric part S' containing the terms 

S cS = S[n} + S'[n,h] . (4.11) 

Due to the invariance of S under global 0(3) rotations it is an 0(3)-singlet and hence 
annihilated by the total angular momentum, 

LS = 0, L = Y,Lx, (4.12) 

X 

such that LS e ff = LS' . Thus, summing over all lattice sites x in (£0]) yields the (broken) 
Ward identity, 

(LF[n] - F[n}LS'[n,h}) =0 , (4.13) 

where the second terms contains the infinitesimal change of the non-invariant part S' of 
the effective action under rotations of n. Note that the coupling constants Xj of the 
0(3)~symmetric operators Sj have disappeared in the Ward identity ( 4.13j ) so that only 



the symmetry-breaking couplings A' are present. We have collected the explicit lattice 
Schwinger-Dyson and Ward identities used in our simulations in App. [C]. As the former 
are local relations, they naturally contain more information than the global Ward identities. 
In particular, one does have access to length scales. 
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5. Comparing Yang— Mills and FN ensembles 



5.1 Leading order ansatz 

To leading order (LO) in the derivative expansion we have a standard nonlinear sigma 
model with symmetry-breaking term, 

S e s = ^2{Xn x ■ An x + X'n x ■ h) , h = e z . (5.1) 

X 

Inverse Monte Carlo amounts to determining the couplings A and A' such that the probabil- 
ity distribution associated with the LO action (5.1) fits the observables of the Yang-Mills 
ensemble of n-fields 2 . The associated Schwinger-Dyson equation ( [4.10] ), with F[n] = n x , 
can be written as 

XH xy + X'G^ y = -Tl8 xy , (5.2) 
where H denotes the (antisymmetrized) two-point function of n % and n'An 3 , 

H xy ee (n*njAnJ) - (4^ An*} ee (n>j/Ar$) . (5.3) 



To analyse (5.2) we define a 'reduced' two-point function h xy and magnetization u, 

h xy = H xy /G^ y , fx = M/G^ x , (5.4) 



and rewrite ( |5.2p as the inhomogeneous system (using translational invariance to replace 

x-y^x), 

\h x + X' = 0, x = l,...,8, (5.5) 
Xh + X' = —fj, . (5.6) 

The solution is found to be 

A= t-^-t-A*, (5-7) 
h x - h 

X' = ~— hjL T~V, (5.8) 
n x — ho 

with the numerical boundary value given by ho = 0.1410 (cf. Fig. ^). Clearly, the system 
J5.5D , (|5.6[) is overdetermined (nine equations for two unknowns). This is reflected in the 
fact that A and A' in (]5.7|) and (|5.8| ) depend on the lattice distance x via h x . If the Yang- 
Mills ensemble were exactly described by the LO action ( |5.1| ), there would be no such 
x-dependence. Rather, for any x = 1, . . . , 8, we would have the same values for A and A', 
respectively. Thus, to test the quality of the LO ansatz, we divide (5^) by (5/7) showing 
explicitly that h x should be constant, 

A' 

h x = — — ee —k' = const. , (5-9) 
A 

Fig. H shows that this is not the case. Therefore, a minimal sigma model with symmetry- 



2 Throughout this section, we refer to algorithm AI. 
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Figure 9: h x for the Yang-Mills ensemble. The boundary value is ho = 0.1410. 



breaking term does not yield a good representation of our Yang-Mills ensemble of n-fields. 
If we nevertheless insist on the LO description, we have to 'fit' h x by a horizontal line 
so that the numerical determination of the couplings via ( |5.7| ) and Q5.8| ) is beset by large 
errors, 



A = -1.41 ± 5.25 , 
A' = -1.33 ±0.74 . 



(5.10) 
(5.11) 



Obviously, A (including its sign) remains essentially undetermined. For A' the situation is 
slightly better. 

In order to assess the errors it is worthwhile to check whether our numerical accuracy 
is sufficient to really validate the Schwinger-Dyson identity (|5.2| ) for the LO action ( |5.1| ) 
on the lattice. To this end we have simulated ( |5.lD with a combination of Metropolis and 
cluster algorithms producing a number of 150 configurations using the central values ( 5.1C| ) 



and (5.11) as the input couplings. The result for h x in the LO ensemble is presented in 



Fig. 10. It is reassuring to note that the simulation of the minimal sigma model reproduces 
the input value «' = 0.943 very well (for i/O), the error being of the order of one percent. 
The prediction (5.9) thus can be verified with high accuracy for the LO action (5.1). We 
conclude that inverse Monte Carlo works quite well when applied to the minimal cr-model. 

The discrepancy between the LO and Yang-Mills ensembles can be further visualized 
by looking at the susceptibility. For the action ( |5.1| ) and the choice F[n] = n", the Ward 
identity ( pL"il|) assumes the simple form 



X 



-SDT/A' 



(5.12) 
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Figure 10: h x for the LO ensemble. Fitting the data points yields h x 



^0.938 ± 0.002. 



A consistency check is provided by noting that this can directly be obtained by summing 
Q5.2| ) over x. Plugging in the magnetization from Table 1 and A' from ( |5,11| ) we find 



X 



0.33 ±0.18 , 



(5.13) 



This is way off the Yang-Mills value of 92.6 displayed in Table 1. For magnetization and 
mass gap the simulation of the LO ensemble yields the values 



SDt = 0.93 , M=1.5, 



(5.14) 



which are both larger than the Yang-Mills values of Table 1. 

The discussion of this subsection thus shows quite clearly that more operators will 
have to be included in order to possibly make inverse Monte Carlo work reasonably well. 

5.2 FN action with symmetry breaking term 

In this subsection we consider the FN action ( |1.1| ) with a LO symmetry-breaking term, 

(n • An) 2 — (n • d^d v n) 2 



S, 



cff 



Xn x ■ An x ± A r 



±A 



'n x ■ /i j 



(5.15) 



This ansatz does not include all terms of next-to-leading order (NLO) in the derivative 
expansion. It should be viewed as a minimal modification of the original FN action by 
adding an explicit symmetry-breaking term to obtain a mass gap. 
The Schwinger-Dyson equation generalizing (f>.2\) becomes 



XH xy ± X FN H™ ± X'G^ 



xy 



-WIS. 



xy ■ 



(5.16) 
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Figure 11: h™ for the Yang-Mills ensemble. The boundary value h Ft 



-31.24 is not displayed. 



The new two-point function H FN is given by (C.19). The local identities (5.16) are to 
be solved for the three unknown couplings A, A' and A FN . Introducing another reduced 
two-point function, 

hxy = Hxy/Gj: y , (5-17) 

which is plotted in Fig. [□], we obtain, instead of ( ^ ) and (|5.6|), the (overdetermined) 
system, 

Xh x + X FN h F x N + A' = , x = l,...,8; (5.18) 

(5.19) 
(5.20) 



Xh y + X FN h FN + X' = , y > x ; 
A/iq + A FN /ig N + A' = -/x , 



to be solved for each pair of lattice distances (x,y), y > x. The number of independent 
pairs is 7(7 + l)/2 = 28 for lattice extension L = 16. The solutions, labelled by x and y, 
are 



A 



'xy 



X — ; 



Afn — 

where we have defined the determinant 



"-xy 
fax hn 



dxy 



, , , = h0(h™ - ^ N ) - hl"(h X - hy) + h X h. 

ihat the numerical values for d xy 
problem with small denominators in the solutions. 



FN 

y 



h y h x 



(5.21) 
(5.22) 
(5.23) 

(5.24) 



We have checked that the numerical values for d xy are not close to zero so that there is no 
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Figure 12: Variation of the couplings ( |5. 21 - ^23 ) with lattice distances x and y. The values 
associated with y > x are plotted along vertical lines labelled by x = 1, . . . , 8. The horizontal lines 
denote the mean value. Note that, for larger distances, the data seem to deviate less from the 
central values, (a) coupling A, (b) coupling A', (c) coupling A FN . 
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For each pair of lattice distances (x, y) we thus have a certain value for any of the three 
couplings ( |5.21 5.23). For each particular coupling those would all agree (within statistical 
errors) if the NLO action ( 5.15| ) would exactly describe the Yang-Mills ensemble. Again, 
however, analogous to the LO case, the couplings do vary with lattice distances x and y as 
shown in Fig. [l^. Numerically, one finds, 



A = -0.232 ± 0.035 , 
A' = 0.257 ± 0.014 , 
A FN = -0.0402 ± 0.0004 . 



(5.25) 
(5.26) 
(5.27) 



Several remarks are in order. First of all, the relative errors, given by the standard deviation 



from the mean (see Fig. 12), are small compared to the LO ansatz. In particular, the signs 
of all couplings are fixed. Interestingly, the addition of the FN coupling A FN , although small 



numerically, has a large effect: it reverts the sign of A' as compared to (5.11), implying a 
negative magnetization. This follows, for instance, from the Ward identity Q5.12j ), which 
still holds for the action ( |5.15 ), and the positivity of the susceptibility, hence 



-x'x ± < o 



(5.5 



in contradistinction with the positive Yang-Mills value of Table 1. 



To further analyse the result for the couplings, we divide ( 5.18 ) by A, leading to a 
linear relation between h and h FN (for i / 0), 



h. 



Afn/A 



(5.29) 



Thus, plotting h x against h™ should yield a straight line with intercept — «' and slope 
— k fn . The numerical values fl5.25 - 5.27 ) yield 



-1.108 ±0.228 , 



0.173 ±0.024 . 



(5.30) 



In analogy with the LO case, we have numerically checked the prediction ( p. 29 ) for the NLO 
action (|5.15| ) by a Monte Carlo simulation with 150 configurations using the input couplings 
( 5,25 - 5^271 ), Fig. [l^ clearly demonstrates the expected linear behavior. A corresponding 
fit results in 

k' = -1.120 , k pn = 0.171 , (5.31) 



being consistent with the central values of ( 5.30| ) to within one percent. We thus conclude 
that inverse Monte Carlo also works quite well for the NLO ensemble. For the sake of 
explicit comparison with Fig.s || and 11 we display the reduced two-point functions ob- 
tained by simulating the NLO action in Fig.s |l4| and [l^. 

As expected, the NLO simulation yields a negative magnetization, 

9Jt =-0.49, (5.32) 

while the mass gap becomes M = 1.2, i.e. slightly larger than the value listed in Table 1. 
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Figure 13: h vs. h FN for the NLO action ( p^ ). 
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Figure 14: /i^ for the NLO ensemble. 



If the Yang-Mills ensemble has anything to do with the FN one, then plotting h 
vs. h™ (as obtained from Yang-Mills) should also show straight-line behavior, at least 
approximately. Fig. 16 displays a linear fit to the the Yang-Mills data with parameters 

k' = -1.232 , k pn = 0.170 . (5.33) 
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Figure 15: h™ for the NLO ensemble. 




Figure 16: h vs. h FN for the Yang-Mills ensemble. The error bars reflect the fact that, for large 
distances x (i.e. x ~ L/2 = 8 for the lattice size used), h and h™ are obtained by dividing two 
small numbers (from the tails of the two-point functions G- 1 , H and H FN ). 



Again, within error bars, these values are consistent with the preceding analysis ( 5.3CQ . 
Note that reverting the sign of A' amounts to reverting the sign of the intercept in Fig. [l6| . 
The data points clearly do not support anything like that. On the contrary, it seems that, 
for small h FN (corresponding to small distances x, see Fig. [TT| ) , the data points deviate from 
a straight line. Playing around with different fits indicates that h rises with a negative power 
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Figure 17: Variation of n' with lattice distances x and y. The central value is 



1.38. 



of h FN for small h FN so that in reality there may be no intercept at all. 

The behavior of function of lattice distances x and y may also be investigated. 

Dividing ( ggp b y (EUD we nnd 



A'/A 



hyh™ 



hi 



h FN 

fly 



-1.38 ±0.29 , 



(5.34) 



where, in the analytic expression, the determinant d xy has dropped out. Fig. |17] shows the 
variation of k' with x and y. Again, a different sign for «' is completely out of reach. 

Following the logic of gradient expansions, one may argue that the effective action 



(5.15) is supposed to represent the Yang-Mills ensemble only for large distances. Fig. 16 
for instance, seems to indicate that the straight-line fit works particularly well for the last 
three points to the right which correspond to x = 6, 7, 8, respectively. In physical units, 
this amounts to distances R larger than six lattice units, i.e. R ~ 0.8 fm. Restricting to 
the analogous data points in Fig. 12, we obtain for the couplings in ( |5.15| ), 

A = -0.2775 ± 0.026 , (5.35) 
A' = 0.2661 ±0.012 , (5.36) 
A FN = -0.040 ± 0.0005 , (5.37) 



and for the 'reduced' ones, 



-0.97 ±0.13 , 



0.15 ±0.01 



(5.38) 



All these do not differ significantly from the values ( 5.25 - 5T27| ) and ( |5.30 ) obtained by 
using the unrestricted data set. In particular, the sign of A' remains positive. We therefore 
conclude that, also at large distances, the minimally modified FN action ( |5.15| ) fails to 
describe the Yang-Mills ensemble of n-fields. 
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6. Summary and discussion 

We have performed a lattice test of the FN conjecture stating that low-energy SU (2) Yang- 
Mills theory is equivalent to a Skyrme-type sigma model. More specifically, FN suggest 
that the knot solitons of their model might be related to the Yang-Mills glueball spectrum. 

Using standard Monte Carlo techniques, we have generated an ensemble of SU (2) 
link fields from the Wilson action. This ensemble was then used to extract an associated 
ensemble of color vectors n, with n parametrizing the gauge invariant distance between 
the maximally- Abelian and lattice Landau gauge slices. As these gauges are close to each 
other, there is a preferred direction for the n-field which corresponds to explicit symmetry 
breaking. In this way we avoid the appearance of massless Goldstone bosons and thus 
generate a nonvanishing mass gap. A study of the exponential decay of correlators yields a 
mass gap close to 1 GeV. To find the effective action describing the Yang-Mills ensemble of 
ra-fields we have employed inverse Monte Carlo techniques. These are based on Schwinger- 
Dyson and Ward identities which we have derived analytically on the lattice. The identities 
have been evaluated numerically for the Yang-Mills ensemble on the one hand, and for 
ensembles stemming from LO and NLO effective actions on the other hand. As a result, 
we have found strong evidence that the ensemble generated from Yang-Mills theory cannot 
be described by the FN action plus a minimal symmetry-breaking term to allow for a mass 
gap. This follows from a number of discrepancies between the two ensembles. First, and 
most prominent, the sign of A' is positive, implying negative magnetization 9Jt, at variance 
with the value from the Yang-Mills ensemble. Second, the reduced two-point function 
h (/i FN ) from the NLO ensemble increases (decreases) with lattice disctance x, while for 
Yang-Mills the behavior is just the opposite. Third, the size of the mass gap is larger than 
for the Yang-Mills ensemble of n-fields. 

It is quite conceivable that magnetization (and susceptibility) can be recovered cor- 
rectly by adding more (symmetry-breaking) terms to the NLO action (work in this direction 
is under way). The same remark applies to the mass gap. Note, however, that one cannot 
naturally expect the Yang-Mills and cr-model mass gaps to coincide due to the nonlocal 
relation between n and the link variables U 3 . Whether this represents a problem is a 
question of scales. If the effective cr-model were valid only for distances of, say, R ~ 0.8 
fm corresponding to energies E ~ 0.25 GeV, as suggested by the discussion of Section 5, 
then it would make no sense to address questions like the glueball spectrum. An analogous 
situation holds for the Fermi theory of weak interactions which also is only effective much 
below the W and Z scales. 

Finally, one should mention that there is still another fundamental problem associated 
with describing Yang-Mills theory in terms of effective cr-models. Allowing for finite tem- 
perature, the latter are in the universality class of the 4d Heisenberg model, while SU(2) 
Yang-Mills theory is in the 3d Ising class This issue has been discussed recently in the 

As there is no unique or natural definition for n, one may try alternative prescriptions for n = n[U]. 
A fairly local one is the following. Write the (gauge fixed) links as U x ,n = u° „ + iUx^. Then define 
n x = u x /\u x \ with the link average u x = X^ Ua; ,n. Under global gauge transformations this transforms 
properly such that n is another color unit vector. 
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context of constructing effective actions via Abelian projections (2^, |3(| [[!). Again, if the 
(T-model scale were below the critical temperature, the effective theory would only be valid 
in the confined phase and would have nothing to say about the behavior of Yang-Mills 
theory close to the phase transition. Otherwise, an infinite number of operators would be 
required which, of course, is anything else but an 'effective' description. Summarizing, we 
conclude that, while a reasonable effective model generalizing the FN action may exist in 
principle, it will be of little practical use. 
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A. Conventions 

Left and right lattice derivatives are defined as 



The ordinary lattice Laplacian A 
on lattice functions / is given by 





(A.3) 



The covariant Laplacian A[C/] in the adjoint representation acts as 




(A.4) 



where we have defined the adjoint link 




(A.5) 



B. Relating LLG and MAG 



From (2.5) it follows immediately that the LLG minimizes the functional p^] 




(B.l) 
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and thus tends to bring the links Ui close to 1. The MAG, on the other hand, minimizes 



Fmag^EC 1 -^ 33 ), 



(B.2) 



and thus wants to bring the 33-entry of the adjoint link Rf close to 1. From (A. 5) it is 
obvious that, if XJ\ equals unity, the same will be true for R\. This can be made more precise: 
if Ui ~ 1 + iaAi, A[ hermitean, then it is an easy exercise to show that Rf 3 = 1 + 0(a 2 ). 
In this sense, the LLG is close to the MAG. 

C. Schwinger— Dyson equations and Ward identities 

We begin with computing the infinitesimal rotations of the various contributions in 
and (^^) to the effective action. It turns out that, for all Sj, the action of the angular 
momentum can be written as 

iL x Sj 



n x x K j x , 



(and analogous for the Si.) with the vectors Kj x and K' kx given by 



(C.l) 



K\ x = 2An x 
K 2x = 2 A 2 n x 

K 3x = 2 [An x (n x ■ An x ) + A (n x (n x ■ An,))] 

dld u n x (n x ■ d ] d u n x ) + d\d u (n x (n x ■ d\d v n x 



Kax — 2 

K' lx = h 

K' 2x = 2h (n x ■ h) 

K' 3x = h (n x ■ An,) + An^n, • h) + A(n x (n x ■ h)) . 



(C.2) 
(C.3) 
(C.4) 
(C.5) 
(C.6) 
(C.7) 
(C.8) 



Choosing the F's in ( 4.13 ) as n x , n % n y and n x n y n c ZJ respectively, results in the Ward 
identities 



^xy A l zu xxy A 2 ^xxy A 3 

^xyz A l ZKj xxyz A 2 ' u ' xxyz A 3 

GU33 \/ I ns~i3ii33 \' i ri&.H33 \l 
xyzz lA \ z(j xxyzz' A 2 ~>~ ^ xxyzz' A 3 



~2G y 



-i33 
'yz 

-(333 



yz ' yz 



(C.9) 
(CIO) 

(C.ll) 



Here, the superscript (i3) denotes symmetrization in i, 3, and we have introduced the 
shorthand notations 



G 



abc. 
xyz.. 

^xy.. 



(n> b y n C z 

E/~iii... 
u xy.. 



(C.12) 
(C.13) 



G 



3ii... 
xxy.. 



^(nlnin y ...) 



X 



G^;^(«An»;...) 



(C.14) 
(C.15) 
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with x denoting summation over all lattice sites x. In particular, one has 



G 



G 



2 X A 



xy — ^^xy 



^xy 



G 



■>-y 



(C.16) 



It seems that a particular recurrence pattern arises in (C.G-G.ll) that could be used to de- 
rive a Ward identity for an insertion of an arbitrary number of n's. For the NLO derivative 
expansion, however, the three identities are sufficient to determine the symmetry-breaking 
couplings A'. 

For the particular case that F in (gTo|) equals n, F£ = n", we can give the general 
Schwinger-Dyson equation in closed form, 



7, ^jHj,xy + ^'kH'k 



xy 



-MS. 



xy , 



(C.17) 



j k 

where we have introduced the two-point function 



H 



3,%y 



(C.18) 



and, completely analogous, H' k . The FN terms are ( |C.4|) minus (|C.5| ), so that the relevant 
two-point function becomes 

(C.19) 



ttFN 
12 xy 



H3,xy H^ X y 



which has been used in ( p. 16 ). 
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